Overlap Equivalence in the Edwards-Anderson model 
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We study the relative fluctuations of the link overlap and the square standard overlap in the 
three dimensional Gaussian Edwards- Anderson model with zero external fleld. We flrst analyze the 
correlation coeflicient and find that the two quantities are uncorrelated above the critical temper- 
ature. Below the critical temperature we find that the link overlap has vanishing fluctuations for 
fixed values of the square standard overlap and large volumes. Our data show that the conditional 
variance scales to zero in the thermodynamic limit. This implies that, if one of the two random 
variables tends to a trivial one (i.e. delta-like distributed), then also the other does and, by conse- 
quence, the TNT picture should be dismissed. We identify the functional relation among the two 
variables using the method of the least squares which turns out to be a monotonically increasing 
function. Our results show that the two overlaps are completely equivalent in the description of the 
low temperature phase of the Edwards- Anderson model. 



The low temperature phase of short-range spin-glasses 
is among the most unsettled problems in condensed mat- 
ter physics 0, • To detect its nature it was originally 
proposed an order parameter by Edwards and Anderson 
[Sj, the disorder average of the local squared magnetiza- 
tion 



qea = Av(a',^) = Av 



(1) 



which coincides with the quenched expectation of the lo- 
cal standard overlap of two spin configurations drawn ac- 
cording to two copies of the equilibrium state carrying 
identical disorder 



Av(Wj^) =< qi >= Av 



-!3{H^+Hr 



(2) 



The previous parameter should reveal the presence of 
frozen spins in random directions at low temperatures. 
While that choice of the local observable is quite natu- 
ral, it is far to be unique; one can consider, for instance, 
the two point function Av(a;,^j ). In the case of nearest 
neighbour correlation function this yields the quenched 
average of the local link overlap. 

When summed over the whole volume, link overlap and 
standard overlap give rise to a priori different global or- 
der parameters. In the mean field case the two have a 
very simple relation: in the Sherrington-Kirkpatrick (SK) 
model, for instance, it turns out that the link overlap co- 
incides with the square power of the standard overlap up 
to thermodynamically irrelevant terms. But in general, 
especially in the finite dimensional case of nearest neigh- 
bor interaction like the Edwards- Anderson (EA) model, 
the two previous quantities have a different behavior with 
respect to spin flips: when summed over regions the first 
undergoes changes of volume sizes after spin-flips, while 
the second is affected only by surface terms. 



From the mathematical point of view, their role is also 
quite different. The square of the standard overlap repre- 
sents, in fact, the covariance of the Hamiltonian function 
for the SK model, while the link overlap is the covari- 
ance for the EA model. Two different overlap definitions 
are naturally related to two different notions of distance 
among spin configurations. It is an interesting question 
to establish if two distances are equivalent for the equi- 
librium measure in the large volume limit and, if yes, 
to what extent (see for a broad discussion on overlap 
equivalence and its relation with ultrametricity). They 
could in fact be simply equivalent in preserving neigh- 
borhoods (topological equivalence) or they could preserve 
order among distances (metric equivalence) . The a-priori 
different properties of the two overlaps have also been dis- 
cussed in relation to the different pictures (droplet 
mean-field 0], TNT [HIlJ), that have been proposed 
to describe the nature of the low temperature spin-glass 
state. In this respect, the distributions of the two over- 
laps are expected to be delta-like (trivial distribution, 
droplet theory), to have support on a finite interval (non 
trivial distribution, mean-field theory), or to have dif- 
ferent behaviour depending on which overlap is consid- 
ered (trivial link overlap distribution, non trivial stan- 
dard overlap distribution, TNT theory). 

In this paper we consider the EA model in d=3, with 
Gaussian couplings and zero external magnetic field in 
periodic boundary conditions. We study the relative fiuc- 
tuations of the link overlap with respect to the square 
of standard overlap. We use the parallel tempering al- 
gorithm (PT) to investigate lattice sizes from L = 3 to 
L = 12. For every size, we simulate at least 2048 disorder 
realizations. For the larger sizes we used 37 temperature 
values in the range 0.5 <T< 2.3. The choice of the low- 
est temperature is related to the possibility to thcrmalize 
the large sizes, but our results are perfectly compatibile 
with those obtained by Marinari and Parisi at T = 
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TABLE I: Parameters of the simulations: system size, number 
of sweeps used for thermalization, number of sweeps for mea- 
surement of the observables, number of disorder realizations, 
number of temperature values allowed in the PT procedure, 
temperature increment, minimum and maximum temperature 
values. 



(see the last paragraph before the conclusions). The 
thermalization in the PT procedure is tested by check- 
ing the symmetry of the probabihty distribution for the 
standard overlap q under the transformation q — » —q. 
Moreover, for the Gaussian couphng case it is available 
another thermahzation test: the internal energy can be 
calculated both as the temporal mean of the Hamiltonian 
and - by exploiting integration by parts - as expectation 
of a simple function of the link overlap We checked 
that with our thermalization steps both measurements 
converge to the same value. All the parameters used in 
the simulations are reported in Tab^ 

We recall the basic definitions. For a 3-dimensional 
lattice A of volume N — L^, the square of the standard 
overlap among two spin configurations cr, t G {+1,-1}^ 
is 



^^^^ 



(3) 



The link overlap is instead obtained from the nearest 
neighbor spins, namely for b — with i,j S A, 

|z — j'l = 1 and (Tb = cridj 



(JbU- 



(4) 



First, we investigate the behavior of the correlation coef- 
ficient between and Q 



P = 



<{q'- <q^>){Q-<Q>)> 
v/< {q^- < q^ >)2 X (g- < Q >)2 > ■ 



(5) 



This quantity will tell us in which range of temperatures 
the two random variables are correlated. In that range 
we further investigate the nature of the mutual correla- 
tion by studying their joint distribution and, in partic- 
ular, the conditional distribution P{Q\q^) of Q at fixed 
values of q^ . We are interested in understanding if a func- 
tional relation among the two quantities exists, i.e. if the 
variance of the conditional distribution shrinks to zero 
at large volumes and around what curve the conditional 



distribution is peaked. We have: 



P{Q\q') = 



Av 



^ CT,X 



AV^ i:„,,e-m«.+« 



(6) 



For this conditional distribution one could compute the 
generic fc— th moment 

Gk{q') := E (g'^lg^) = j'^ Q'^P{Q\q')dQ. (7) 



We will be interested in the mean 

G{q'):^G,{q')=E{Q\q^) 
and the variance 

Var{Q\q')=G2{q')-Gl{q'). 



(8) 



(9) 



The method of the least squares immediately entails that 
the mean G{q^) is the best estimator for the functional 
dependence of Q in terms of q^. In fact, given any func- 
tion ^.(g^), the mean of {Q — h{q'^))'^ according to the joint 
distribution P{Q,q^) is " h{q^)fPiQ,,q^) = 

Ej P{q]) E^iQ^ - Hq])?pIQMj), where the sums run 
over all possible values of the random vector {Q,q^)j 
which are finitely many on the finite system we simulated. 
Therefore, to minimize the mean it suffices to minimize 
the inner sum, i.e. to choose h{q^) as the mean G{q^) of 
Q with respect to the conditional distribution 

The scaling properties of the conditional variance © 
and the functional dependence (O provide important 
informations about the low temperature phase of the 
model. Indeed, a vanishing variance in the thermody- 
namic limit implies that the two random variables Q and 
q^ do not fiuctuate with respect to each other. If the 
functional dependence G{q^) among the two is a one-to- 
one increasing function, then it follows that the marginal 
probability distributions for the standard and link over- 
lap must have similar properties. In particular, if one 
of the two is supported over a point then also the other 
must be so. 

We now describe our results. Fig. ^ shows the corre- 
lation between the square standard overlap and the hnk 
overlap. The plot of Eq. is done for different sizes 
of the system as a function of the temperature. It is 
clear from the figure that, as the system size increases, 
the correlation remains strong in the low-temperature re- 
gion, while it becomes weaker in the high temperature 
region. A sudden change in the infinite volume behavior 
of p can be expected to occur close to the critical temper- 
ature Tc of the model. The best estimate available in the 



3 



literature - obtained through the analysis of the Binder 
parameter's curves of the variable for different system 
sizes - gives = 0.95 ±0.04 01 (we independently repro- 
duced this estimate with our data for and obtained an 
estimate of = 0.95 ± 0.03). 

For each temperature, we did a fit of the data for p to 
the infinite volume limit. We tried different scaling for 
the data, both exponential pl{T) = Poo(T) + a(r)e''(^)^ 
and power law pl{T) = /9oo(T) + a{T)L^('^\ The inter- 
esting information is contained in the asymptotic value 
Poo{T). We measured the normalized for different 
values of Poo{T) in the range [0, mini, pi(T)] and keep- 
ing a{T) and h{T) (or a{T) and (i{T)) as free parame- 
ters. In the region T > 1.0, we found that attains 
his minimal value for Poc{T) = 0. For T < 0.9 the 
develops a sharp minimum corresponding to values 
Poo{T) ^ 0. The whole plot of the curve Poo{T) as ob- 
tained from the best fit is represented in Fig. Also, 
in the inset of the same figure it is shown the standard 
finite size scaling of the data. We plot pl{T)L'^^'' versus 
the scaling variable {T ~ Tc)L^^'^ . A good scaling plot 
is obtained using = 0.95, v = 0.71, ?A = 0.038. The 
discrepancy between the value v — 2.0 of ref. ^ has to 
be attributed to the non-linear relation between Q and 
(see below). Fig. Q tells us that in the high tem- 
perature phase the two random variables standard and 
link overlap are asymptotically uncorrelated while in the 
low temperature one they display a non- vanishing corre- 
lation: within our available discrete set of temperature 
values, the temperature at which the correlation coeffi- 
cient starts to be different from zero is in good agreement 
with the estimated critical value of the model. 

We consider then the problem of studying the func- 
tional dependence (if any) between the two random vari- 
ables Q and q^ in the low temperature region. The points 
in the Fig. Q show the function G(g^) of Eq.® for dif- 
ferent system sizes at T = 0.5, well below the critical 
temperature. Also we studied a third order approxima- 
tion of the form Q — g[q^) ~ a + bq^ + cq^ + dq^ . Since 
we must have Q — \ for = 1, this actually implies 
d = 1 — a — b — c. The coefficients aL,T, bL,T, cl,t have 
been obtained by the least square method and then fit- 
ted to the infinite volume limit. The result is shown 
as continuous fines in Fig Q. The good superposition 
of the curves to the data for G{(f') indicates that the 
functional dependence between the two overlaps is well 
approximated already at the third order. 

Finally, we measured the variance Eq. © at low tem- 
peratures. We observed that the distribution is concen- 
trating for large volumes around its mean value. The 
trend toward a vanishing variance for infinite system sizes 
is very clear. We analyzed all temperatures below Tc and 
we found that the best fit of VarL{Q\q^), in terms of the 
X^, is obtained by a power law of the form a{T)L^''^'^'> +c 
and it gives c = for every value of the temperature. For 
the lowest available temperature T = 0.5, this is shown 



in Fig |@J where we plot the data for VarLiQ\q^)L^'^^ 
for different system sizes L: all the different curves col- 
lapse to a single one. The data for other temperature 
values behave similarly, the only difference being that 
the coefficient b{T) is increasing with the temperature 
(it stays in the range [1.43, 1.74] for T e [0.5,0.9]). This 
result has quite strong consequences because it says that 
the two random variables Q and q^ cannot have different 
triviality properties: if one of them is trivial (delta-like 
distributed) the scaling law for their conditional variance 
implies that also the other is trivial. Our result rules then 
out the possibility to have a non-trivial standard overlap 
with a trivial link overlap as predicted for instance in the 
so called TNT picture U^^M 

It is interesting to compare our result with previous 
works. Marinari and Parisi have studied the relation 
Q = (1 - A{L)) + {A{L) - B{L))q^ + B{L)q'^ among the 
two overlaps at zero temperature, by ground state per- 
turbation. We have extrapolated our data in the low 
temperature regime to zero temperature by a polyno- 
mial fit and then to the infinite volume limit {L = oo). 
The best fit for L = oo (i.e. the one with smaller x^) is 
quadratic in L'^. It gives A = 0.30 ± 0.05 (x^ = 0.21), 
which is in agreement with the independent measure of 
Marinari and Parisi {A = 0.30 ± 0.01, x^ = 0.6). Note 
that their results are obtained with a complete differ- 
ent method than Montecarlo simulations, namely exact 
ground states. Sourlas studied the same problem in a 
different setting called soft constraint model. Although 
a direct quantitative comparison is not possible with our 
method, his results are qualitatively similar. In the con- 
text of out-of-equilibrium dynamics, a strong correlation 
between link-overlap and standard-overlap in the low- 
temperature phase was pointed out in ref. 

In conclusion, our result shows quite clearly that, 
within the tested system sizes, the study of the two 
order parameters, the square of the standard overlap and 
the link overlap, are equivalent as far as the quenched 
equilibrium state is concerned. In view of our result, the 
proposed pictures which assign different behaviour to the 
two overlap distributions, in particular the TNT 0iQ] 
picture, should be rejected. It is interesting to point 
out that, since the present analysis deals only with the 
distribution of P{q^,Q) and not with the higher order 
ones like for instance P(gi 21 3i Qi.Zi ^2,3), our results 
are compatible with different factorization properties of 
the two overlaps like those illustrated in [l3| . 

Acknowledgments. We thank Enzo Marinari, Gior- 
gio Parisi and Federico Ricci Tersenghi for interesting 
discussions. We thank a referee for the suggestion to im- 
prove the Fig. (O by a finite scale analysis. We also 
thank Cineca for the grant in computer time. C. Gia- 
rdina acknowledges nwo- project 613000435 for financial 
support. 



4 



[1] M. Mezard, G. Parisi, M.A. Virasoro, Spin Glass Theory 

and BeyondWorld Scientific, Singapore (1987). 
[2] CM. Newman and D.L. Stein, cond-mat/0503345 
[3] S.F. Edwards and P.W.Anderson, Jou. Phys. F., 5, 965, 
(1975). 

[4] G. Parisi, F. Ricci-Tersenghi, J. Phys. A: Math. Gen. 33, 
113, (2000). 

[5] P. Contucci, J. Phys. A: Math. Gen. 36, 10961, (2003). 
[6] E. Marinari, G. Parisi, Phys. Rev. Lett. 86, 3887, (2001). 
[7] E. Marinari, G. Parisi, and J.J. Ruiz-Lorenzo, Phys. Rev. 

B 58, 14852 (1998) 
[8] S. Jimenez, V. Martin-Mayor, G. Parisi, A. Tarancon J. 

Phys. A: Math, and Gen. 36, 10755 (2003) 
[9] N. Sourlas, Phys. Rev. Lett. 94, 70601, (2005). 
[10] F. Krzakala and O. C. Martin, Phys. Rev. Lett. 85, 3013 

(2000) 

[11] M.Palassini, A.P.Young, Phys, Rev. Lett. 85, 3017, 
(2000). 

[12] D.S. Fisher and D.A. Huse, Phys. Rev. Lett. 56, 1601 
(1986) 

[13] P.Contucci, C.Giardina, Phys. Rev. B 72, 14456, (2005). 



*f « g . 



* L=3 
L=4 
L=5 



L=IO 
L=I2 



0.4 0.6 0.8 1 1,2 14 1.6 1 



2 2.2 24 



FIG. 1: phiT) as a function of the temperature T 
for different sizes L of the system. 
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FIG. 2: poo(r) as a function of the temperature 
T. In the inset it is shown the finite size scaling. 
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FIG. 3: Plot of the curves g{q^) (continuous lines) 
and of G{q^) (dotted lines) together with the infi- 
nite volume limit curve gac{q^) (upper continuous 
line) for T = 0.5. 




FIG. 4: L^-'^^VarL{Q\q^) as a function of for 
the temperature T = 0.5 and for difi^crent sizes L. 



